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ABSTRACT 

The study of the dynamical evolution of massive black hole pairs in mergers is 
crucial in the context of a hierarchical galaxy formation scenario. The timescales 
for the formation and the coalescence of black hole binaries are still poorly con- 
strained, resulting in large uncertainties in the expected rate of massive black 
hole binaries detectable in the electromagnetic and gravitational wave spectra. 
Here we review the current theoretical understanding of the black hole pairing 
in galaxy mergers, with a particular attention to recent developments and open 
issues. We conclude with a review of the expected observational signatures of 
massive binaries, and of the candidates discussed in literature to date. 



Introduction 



Understanding the formation and evolution of massive black holes (MBHs) is one 
of the most exciting goals of contemporary astrophysics and co smology. It is no\y well 



Ferrarese fc Ford 



20051) 



established that MBHs are ubiquitous in nearby spheroids (e.g. 
most of them lurking in a quiescent accretion state, while during the cosmic history 
million-to-billion solar mass MBHs powered quasars. These objects have become in the 
last years central building blo cks for all the proposed scenarios of galaxy formation (e.g. 



Kauffmann fc Haehnelt 



20001). playing a major ro le in shaping galaxies through feedback 



processes (jSpringel. Di Matteo fc Hernquist 



formed during the galaxy merging process ( iBegelman. Blandford fc Rees 



200511. Massive black hole binaries (MBHBs 



19801 ) ■ promise 



to be among the most luminous gravitational wave (GW) sources for future space-borne 
interferometers like the proposed New Gravitational wave Observatory (NGO) Q, and 



^https: / /lisa-light . aei.mpg.de /bin / view / 
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ongoing Pulsar Timing Array (PTA) campaigns ( iHobbs et al. Il2010l ). Theoretical modeling 
of MBHB dynamics is essential in addressing a number of fundamental astrophysical 
questions (such as the merger-quasar connection or the MBH-host relations), and in 
identifying putative signatures that my serve as a guidance for present and future 
observational campaigns. 

Early stages of MBH pairing have been observed, from the initial phases of galaxy 



mergers, where two distinct bu ; 



of ~ 100 kr 


1 

)c (IHennawi et al. 


2006; 


Shen et al. 


2010 




Farina et al. 


2011) 



gravi t ationally bourid ga. 



Myers et al. 



2007 



axies are observable at separa tions 



2008 



Foreman et al. 



2009 



20111 ). down t o unbound pairs of MBHs at separations o 



;^ 1 kpc embedded in a single galaxy remnant (IKomossa et al. 



2003 



Colpi fc Dotti 



201l|). 



During this initial stage the MBH pairing is driven by dynamical friction acting on the 
host galaxies. The two MBHs (hereafter Mi and M2 for the primary and secondary MBH, 
respectively) bind in a binary if they reach a relative separation 



obhb 



2a2 



0.2 Mbhb,6 o-ioo PC, 



1 



where a is the velocity dispersion of the host galaxy, Mbhb = Mi + M2 is the total mass 
of the binary, and aioo and Mbhb,6 are in units of 100 km s~^ and 10® Mq, respectively. 



a relation (see, e.g. 



3, the mass of t 


le MB 


3B scales with a fol 


Gultekin et al. 


2009 


I 


Graham et al. 


2011 



owing the MBH mass vs. 
, and references therein). 



equation [1] implies obhb ~ 0.5 MgVrg g pc 



The efficiency of the process depends on the mass ratio between the merging galaxies. 
While equal mass (major) mergers result in a fast formation of a MBHB, in very unequal 
mass (minor) mergers the satellite, tidally disrupted along the course of the encounter, 
leaves its naked MBH wandering in the outskirts of the most massive galaxy (see, e.g. 



Governato. Colpi fc Maraschi 



I994J ) . Recent numerical studies have addressed the efficiency 



of the formation of a MBHB as a function of the mass ratio between the two galaxies, the 
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gas fraction in the merging galaxies, and the redshift of the merger (jCallegari et al. Il2009[ ) 
In particular, the presence of a significant gas component in the satellite helps the MBHB 
formation: during the first pericenters, the interaction between the two galaxies promotes 
the formation of bars, that convey a large fraction of the available gas in the center of 



the merging galaxies 


(already noticed in 


Barnes & Hernquist 


1996 




Barnes 


2002) 



2002I ). This new nuclear gas overden sity deepens the 



poten tial well of the secondary nucleus and prevents its tidal disruption. 



Callegari et al. 



( 120091 ) found a critical galaxy mass ratio for the formation of a MBHB of ~ 1/10 for gas 



rich galaxies and > 1/4 for gas poo r galaxies. In zero-order approximation, assu ming the 



MBH mass vs. bulge mass relation ( iHaring fc Rix 



2004 : 



Marconi fc Hunt 



20031 ) we expect 



similar mass ratios q = M2/M1 (< 1) in MBHBs. However, q can be increased for gas rich 
galaxies. The strong gas inflows in the satellite (more perturbed) galaxy can result in a 



faster gr owth of the sma. 



(g > 1/3 



Callegari et al. 



er M BH, and in higher values of the expected MBH mass ratio 
2OI1U!. 



As a consequence, efficient dynamical friction promotes the formation of MBH binaries 
with similar mass ratios. The expected number of observable binaries, and the rate of MBH 
coalescences, however, depends on the dynamical evolution after the binary formation. In 
order to coalesce through GW emission in less than an Hubble time, the two MBHs have to 
reach a separation 

^1/4 / Ml 



acw 2 X 10-7(e 



'BHB 



3/4 



pc. 



(2) 



(l + g)i/2 ^1O6M0 

where /(e) = [1 + (73/24)e^ + (37/96)e^](l — e^)^'^/^ is a function of the binary eccentricity 



^We caution that these ranges in q implicitly assume similar morphology in the two 
merging galaxies, although minor mergers can involve very different morphological types. A 
study of MBH pairing and the formation of MBHBs in mixed mergers (e.g. disks merging 
with ellipticals) is not available to date. 
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e (jPeters Ill964l ). Circular equal mass binaries can coalesce shrinking by a factor 



/ T^/r \ -1/4 



while this factor decreases to 1 for 10^ M© binary with e ~ 0.999. In order to understand 
the final fate of a MBHB and to constrain theoretically its observability, it is fundamental to 
study at the same time the orbital decay and the eccentricity evolution of the binary. The 
final fate of the MBHs, i.e. if they will coalesce in a single object or not, strongly depends 
on the amount of matter (stars and gas) they can interact with after the binary formation. 
A definite answer is not present to date. In this paper different scenarios will be discussed, 
depending on the nature of the galaxy mergers and on the properties of the galaxy nuclei. 

This paper is organized as follow: In Section [2] we review the dynamical evolution 
of MBHB in gas poor environments, while the effect of gas is discussed in Section [31 In 
Section m we describe the MBHB candidates observed to date. Finally, our conclusions are 
drawn in Section [51 



2. Dynamical evolution in gas-poor environment 

In systems where dynamical friction is efficient in dragging the two MBHs to the 
center of the merger remnant, the now bound MBHB is inevitably embedded in a gas and 
star rich environment. Such rich ambient provides a variety of physical mechanisms to 
efficiently extract the energy and angular momentum of the MBHB, promoting its final 
coalescence. In this section, we focus on dynamical processes involving interactions with 
stars. MBHB s in pure stellar environments were the first to be examined (the basics 



going back to 



Begelman. Blandford fc Rees 



1980l ). for the obvious reason that stars can 



be considered as point particles, affected by gravitational forces only. The MBHB-star 
interactions are therefore adequately described by Newton's laws only, without all the 
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complications involved in gas dynamics. Nonetheless, a single star-binary interaction is, 
by definition, a three body problem, and the dynamics of the system is inevitably chaotic. 
Therefore, no simple analytical solutions are viable, and numerical studies (both involving 
three body scatterings and full N-body simulations) have been massively exploited to tackle 
the problem. The fate of the MBHB is determined by its semimajor axis and eccentricity 
evolution (see the introduction); in the following we discuss them separately. 



2.1. Shrinking of the binary semimajor axis 



The basic physical process driving the MBHB evolution in presence of stars is the 
slingshot mechanism. A star intersecting the MBHB orbit undergoes a complex three-body 
interaction being eventually ejected at infinity, carrying away en ergy and angular moraentum 



from the binary. Extensive three body scatterin g experiments (IMikkola fc Valtonen 



Quinlan 



1996 



Sesana. Haardt fc Madau 



1992 



20061 ) have shown that ejected stars carry away 



an energy per unit mass of the order of 



3 G/i 



BHB ~ 77- 



(4) 



2 CtBHB 

where /z = M1M2/MBHB is the reduced mass of the binary. Assuming a classical interaction 
rate given by F = p/mSf, where S is the binary cross section, p/m = n is the number 
density of the ambient stars, and v is their typical velocity ' at infinity ' with respect to the 



binary (i.e. the velocity dispersion a of the stellar system), I Quinlan 
the evolution of the binary semimajor axis is simply given as 



(119961 ) showed that 



da 



BHB 



'-BHB 



Gp 



dt 



a 



(5) 



where if is a dimensionless hardening rate. If oehb < GM2/(4(T^), if ^ 15 independently 
on the mass, mass ratio and eccentricity of the system. In principle, given a stellar system 
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with density p and velocity dispersion cr, equation ([5]) predicts efficient coalescence of the 
MBHB. 

However, the above simple treatment ignores the concept of loss cone depletion. In 
an extended stellar system, only a tiny fraction of the stellar phase space allows orbits 
intersecting the MBHB, commonly referred as 'binary loss cone'. As stars are ejected, the 
loss cone is depleted, and the binary evolution is go verned by the rate at which new stars 



are fed into the loss cone (iMakino fc Funato 



20041 ). In typical stellar systems, the mas s 



in sta rs in the binary loss cone is of the order of few times ^ fjMerritt fc Milosaylievic 



20051 ). insufficient to reach ag^ in most of the c ases (ISesana. Haardt &: Madau 



This is the origin of the 'last parsec problem' ( iMilosavljevic fc Merritt 



20071 ) 



200l[ ). In a 



spherical stellar system, th e loss cone reffiling proceeds on a two body relaxation timescale 



( iBinney fc Tremaine 



19871 ). which is usually much longer than t 



le Hubble time. In the 



last decade, this fact has beeii confirmed in N-body simulation s (IMilosavljevic fc Merritt 



2001 



Makino fc Funato 



2004 



Merritt. Mikkola fc Szell 



20071 ). In such simulations, after 



loss cone depletion, further hardening was provided by two body relaxation. This is a 
process that depends on the 'granularity' of the systems, and the result is an N-dependent 
hardening rate, with the binary evolution slowing down as the number of particle in the 
simulation increases. Extrapolating these results to a realistic N representative of a galactic 
bulge, the binary evolution would have stalled. 

In recent years, evidence has emerged that the 'last parsec problem' might be 
an artificial product of the 'spherical cow' approximation which is often exploited in 
astronomy. Basically, the spherical systems studied in the simulations represent a worse 
case (and unrealistic) scenario. MBHBs are infact produced in galaxy mergers, in which 
the resulting stellar bulge is rotating, triaxial and likely to undergo bar-like instabilities. 
In a triaxial potential, an orbiting star does not conserve any of its angular momentum 
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components (IBinney fc Tremaine 



19871 ). As a result, there is a vast fan aily of orbits 



(called centrophilic) that are allowed to get arbitrarily close to the binary (jPoon fc Merritt 



2004 



Merritt fc Poon 



2004 



Merritt &: Vasiliev 



201ll ). keeping the loss cone full during 



the MBHB ha r denin g process. Recent N-body simulations have confirmed this scenario. 



Berczik et al. 



( 120061 ) studied the evolution of a MBHB in a rotating bulge. In this case, 
the stellar system experiences a bar instability resulting in a triaxial potential. The 
binary hardening rate was found to be N-independent; a proof that the hardening was not 
proceeding because of spurious two body relaxation. More recently, the advent of GPU 
computing made possible to simulate 'ab initio' the evolution of two interacting stellar 
bulges hosting MBHs; a first s tep toward a realistic galactic merger scenar io. Several 



simulations were perfo r med 



Gualandris fc Merritt 



by 



Khan. Just fc Merritt 



mm 



Preto et al. 



(1201 ih and 



(120111 ). In all cases, the stellar remnant was triaxial and rotating, 
and the hardening rate, given by triaxiality driven loss cone replenishment was found 
to be independent on N, implying coalescence timescales of ~ lO^yr. Remarkably, when 
normalized to the merging galaxy properties, t he binary hardening rates fou nd in these 



simulations follow equation where H ^ 20 (iGualandris fc Merritt 



20111 ). This is a 



consequence of the fact that, whatever is the geometry of the system, the average 'quantum' 
of energy taken away from an interacting star is always the same, and the evolution of the 
system is determined by the star-binary interaction rate only. 



2.2. Eccentricity evolution 

As pointed out in the introduction, eccentricity plays an important role in driving 
the binary coalescence. However, addressing the eccentricity evolution of the system is 
more complicated because the Acbhe caused by each individual interaction depends on 
a combination of energy and angular momentum exchanges. The angular momentum 



- 9 - 



distribution of the interacting stars is therefore crucial. The eccentricity evolution can be 
described as 

K (6) 



BHB 



(iln(l/aBHB) 



Here 7^ is a dimensionless p arameter that, 



ratio, and eccentricity itself (IQuinlan 



differently than H, depends on th e binary mass 



1996 



Sesana. Haardt fc Madau 1120061 ). In general K 



is a positive number in the range — 0.3 (the peak value occou rs at eRHR ~ .6), meaning 



Sesana 



( I2OIOI ) constructed 



that the binary eccentricity grows during the shrinking process, 
complete binary evolutionary tracks by coupling three body scattering experiments of 
bound and unbound stars to an analytical description of the stellar distribution and on the 
loss cone refilling. As a general trend, quasi circular-equal mass MBHBs experience just 
a mild eccentricity growth, while systems which are already eccentric at the moment of 
pairing, or with q significantly lower than 1, can evolve up to cbhb > 0.9. 

The eccentricity evolution in stellar environments has been tackled by several 
authors by means of full N-body simulations. However, the limited number of particles 
(A^ < 10^) in such simulations results in very noisy behavior for the binary eccentricity, 
and it is difficult to draw conc l usions about the general trends behind the numerical 



noise. 



Milosavljevic fc Merritt 



( 120011 ) carried out numerical integration of equal MBHBs 
embedded in two merging isothermal cusps (p cx r""^, with 7 = 2). Starting with 
circular orbits they find a mil d eccentricity increase to a value of < 0.2 during the 



stellar driven hardening phase. 



Merritt. Mikkola fc Szell 



embedded in Dehnen density profiles (jPehnen 



mm considered equal MBHBs 



19931 ) with 7 = 1.2 with different 



initial eccentricities. Again, they find that circular binaries tend to stay circular, while 
eccentric binaries tend to increase their eccentricities in reasonable agreement with the 
prediction of scattering experiments. Simulation s start ing b efore the f ormat ion of MBHBs 



carried by 



Hemsendorf. Sigurdsson fc Spurzem 



mm and 



Aarseth 



mm produce 



binaries with cq ~ 0.8 at the moment of pairing, with e subsequently increasing up 
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to > 0.95. lAmaro-Seoane. Miller fc Freitag I (|2009[ ) focused on binaries of intermediate 
MBHs (M ~ IO^Mq) in massive star clusters. Coupling full N-body simulations 
to three body scattering experiments, they find binaries with significant eccentricity 
(~ 0.5 — 0.6) at the moment of pai ring, growing up to > .9 du ring the hardening phase; 



similar conclusions are reporte d by 



Amaro-Seoane et al. 



f|2010l). On the small q side. 



simulations were performed by iBauni: 



Matsubayashi. Makino &: Ebisuzaki 



gardt. Gualandris fc Portegies Zwart 



(120061) and 



(120071 ). assuming a stellar density profile 7 = 1.75, 



motivated by t he analytical equi 



massive object (IBahcall fc Wolf 



ibrium solution for a dense relaxed stallar cusp around a 



19761 ). When properly rescaled, the eccentricity increase 
found in both papers agrees rema r kably well with predict ions based on the hybrid model by 



Sesana. Haardt fc Madau 



(I2008h . 



Iwasawa et al. 



(I2OIII ) investigated in detail the angular 



momentum exchanges between the binary and the stars responsible for the eccentricity 
growth, in bound stellar cusps. In particular they showed that stars counterrotating with the 
binary tend to extract a lot of angular momentum from the MBHB, causing the eccentricity 
growth, whereas corotating stars do not. This is a sim ple consequence of angular mo r nentu m 



conservation during the ejection process, as shown by 



Sesana. Gualandris fc Dotti 



(1201 ih 



The evolution of t he binary ec c entric ity can be extremely different for non isotropic 



systems. For example. 



Dotti et al. 



(120071 ) showed that at large scales before the formation 



of a binary dynamical friction exerted by rotationally supp orted stellar disks tend to 



circula rize the orbit of a MBH pair. At smaller separations. 



Sesana. Gualandris fc Dotti 



( I2OIII ) demonstrated that in a rotating stellar system, the eccentricity evolution of unequal 
MBHB is dramatically affected by the level of co/counter rotation of the stellar distribution 
with respect to the binary, with corotating distributions promoting circularization rather 
than eccentricity growth. Nonetheless, most of the simulations inv olving rotating bulges 



( Berczik et al. 



2011 



2006 



Preto et al. 



Berentzen et al. 



2011 



2009h. or mer ging systems ( iKhan. Just fc Merritt 



Gualandris fc Merritt 



2OIII ) find quite eccentric binaries at the 
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moment of pairing (ranging from 0.4 to 0.8), and the subsequent evolution leads to a 
general eccentricity growth, in good agreement of what predicted for an isotropic stellar 
distribution. This may be because the binary evolution is mostly driven by loss cone 
refilling of unbound stars on almost radial orbits, with negligible initial angular momentum. 

Overall, the emerging general picture favors efficient coalescence of MBHBs in dense 
stellar merger remnants. The triaxial and rotating nature of the stellar distribution 
promotes efficient loss cone refilling, while large eccentricities (especially in unequal mass 
systems) shorten the gap between the binary pairing and the efficient GW emission stage. 
In the near future, massive N-body simulations with several million particles will offer a 
unique opportunity to confirm this scenario. 



3. Dynamical evolution in gas-rich environment 



3.1. Formation of a MBH binary in a circumnuclear gas disks 



As discussed in the Introduction, in comparable mass, gas rich galaxy mergers 
the gravitational interaction drives strong inflows of gas toward the galactic centers. 



The numerica 



Mayer et al. 



mult i-scale investigation of an equal mass galaxy merger discussed in 



20071 1 revealed that, in advanced stages of the galaxy merger the two MBHs, 
orbiting in the central 100 pc of the merger remnant, are embedded in a dense, rotationally 
supported, gas disk (see flgure[T]). This circumnuclea r disk is self-gravita ting, and can be 



up to ~ 500 times more massive than the MBH pair (IMayer et al 



20071 ). The dynamical 



evolution of the two MBHs is driven by dynamical friction, and, since the circumnuclear 
disk is the densest structure in the remnant nucleus, it is the main cause of their orbital 



decay. 



Mayer et al. 



(120071 1 followed the evolution of the two MBHs from the initial stages 



of the galaxy merger down to to < 5 pc where they form a binary, as the mass of gas 
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enclosed within their separation is less than the mass of the binary. 



The evolution of MBHs in circumnuclear disks has been studied in details in dedicated 
simulations, in which the former evolution of the MBHs at distances ^ 100 pc is not 
explored. This allows to achieve a better resolution in the central region of the remnant 
and to study the la tter MBH pairing . Simi l ar indi pendent inv estigations discussed in 



Escala et al. 



fl2005h and 



Dotti et al. 



(12007 



2009al ) agree with 



Mayer et al. 



(120071 ) on a 



rapid (on a timescale of ^ 10^ yr) formation of a MBHB at parsec separations. The higher 
resolution in these studies allow for a further decay of the binary, that reaches sub-pc 
separations comparable to the spatial resolution in these numerical studies. 



During this intermediate stage (100 pc 



> a > 



0.1 pc) the in teraction with the 



circum nuclear disk strongly affects the eccentricity of the MBH pair. 



Dotti. Colpi fc Haardt 



(120061 ) first noted that, for MBHs with an initial significant eccentricity the decay phase is 
preceded by a circularization of the orbit. This is due to the dynamical friction exerted by 
a rotating background onto the BHs at their apocenters. Since the circumnuclear disk is 
rotationally supported, the BHs at their apocenter move slowlier than the local gas, and are 
dragged in the direction of their motion. This positive torque is exerted only in proximity 
of the apocenter, where the angular momentum of the MBHs can maximally increased. 
This effect is a gen eral feature of dynamical friction exerted by a rotating background. 



Dotti et al. 



(120071 ) demonstrated that the same effect is pr esent for MBHBs orbiting in a 



Callegari et al. 



(120091 ) found 



stellar circumnuclear disk (as discussed in Section [2^21) . and 
the same effect acting in unequal mass galaxy mergers at larger scales, where the disk of 
the primary galaxy circularizes the orbit of the satellite. 
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Mayer et al. 



2003). The 



Fig. 1. — Different stages of an equal mass gas rich merger (from 
color code refers to gas density. The four panels to the left show the large-scale evolution 
at different times. The strong inflow of gas onto the two galactic centers is observable after 
the flrst pericenter. The boxes are 120 kps on a side (top) and 60 kpc on a side (bottom). 
The panels on the upper right are zoom-ins (8 kpc on a side) of the central region during 
the merger of the two galactic cores. The two bottom panels show the circumnuclear disk 
forming in the center of the galaxy remnant. The MBHs are shown as black spheres. 
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3.2. Evolution of close MBH binaries in circumbinary gas disks 



After the formation of a close binary, the MBHs, acting as a source of angular 
momentum, exerts a tidal torque that inhibits the gas from drifting insi de its orbit. This 



creates a hollow density region, called gap, that surrounds t 


le binary (e.g. '. 


_/in & PaDaloizou 


1979a 


b; 


Artvmowicz & Lubow 


1994; 


Sver & Clarke 


1995; 


Gould & Rix 


2000 


)• 



As a consequence of disk clearance, corotation and inner Lindblad resonances 
are reduced in power, drastically changing the dynamical evolution of the binary. 
The same transfer of angular momentum that keeps the disk from infalling onto the 
MBHs drives the shrinking of the binary. Analytical and numerical studies agree 
in predicting that the interaction between the binary and the surrounding material 



reduces the semimajor axis of the binary, and increases the eccentricity of an initially 



low eccentricity binary (e.g. Papaloizou & Prinele 


1977; 


Goldreich & Tremaine 


1980; 


Ivanov. Papaloizou & Polnarev 1999; Gould & Rix 


2000; 


Armitaee & Nataraian 


2002 


Goldreich & Sari 2003; Armitaee & Nataraian 200 


5; Haiman. Kocsis & Menou 


2009; 


Cuadra et al. 2009; Roedie et al. 2011). 



The MBHs/ circumbinary disk interaction can be studied in two limit cases: i) assuming 
that the MBHB is embedded in a virtually infinite disk, as in the case in which it is 
continuously refilled by a long-lived larger-scale structure, or ii) assuming that the disk has 
a finite mass (and, as a consequence, a finite reservoire of energy and angular momentum). 



In the first limit of a constantly fuelled disk, analytical models of the evolution of the 
binary are available. The orbital decay timescale is: 

_ Mdisk(Q) + M2 

f decay , ^ / \ J- vise ; \' , 

Mdisk{a) 

Where Mdisk('3.) is an estimate of the disk mass inside the orbit of the secondary MBH. 
As long as Mdisk(o.) is greater than M2, the MBH behaves as a fiuid element, decaying 
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Fig. 2. — Two consecutive stages of the MBHB decay in the nucleus of a gas rich galaxy 
merger remnant. Left panel: Central region (~ 5 pc) of a large-scale (100 pc) massive 
(lO^M^;^) circumnu clear disk, embedding a sub-pc MBHB (from one of the simulations in 



Dotti et al. 



2009al ). The color code refers to the gas density (high densities in yellow). Right 



pane l: an eccentric MBHB embedded in a sub-pc circumbinary disk (from 
20111). Here darker blue refers to higher densities. 



Roedig et al. 
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on the viscous timescale tvisc onto the primary. As the binary orbit shrinks, Mdisk(a) 
decreases. When the enclosed mass becomes comparable with the secondary its decay 
timescale increases to result ing in very large migration time s cales at small separations. 



However, as first noticed in 



Ivanov. Papaloizou fc Polnarev 



(Il999l ). in a continuously 



refilled disk the stalling of the binary would cause a staedy increase of the density of 
the inner edge of the circumbinary disk, until the mass close to the binary becomes 
again comparable to the second a ry. A t this stage, fast migra tion starts again. Following 



Ivanov. Papaloizou fc Polnarev 



fll999h 



Cuadra et al. 



(120091 ) estimated that, for a disk on 



the verge to undergo fragmentation (i.e. as dense as possible), the coalescence timescale 
of a 3 X IO^Mq binary with q = 1/3, starting from an initial separation of < 0.05 pc, is 
< 10^ yr. This timescale decreases with the decrease of the binary mass. Promisingly, the 
initial separation assumed here is close to the limit achieved in the larger scale simulations 
discussed in Section 13.11 Note, however, that such a timescale is comparable with the age of 
the Universe at z > 7. If the migration in a dense circumbinary disk is the fastest process 
driving the MBHs coalescence, no coalescences of binary with Mbhb ~ lO^M© are expected 
at 2; > 7. 

In a similar way, the s t udy of the evolution of e is po ssible. For simple alpha disks 



jShakura fc Sunyaev 



I973I) 



Artymowicz fc Lubow 



(I1994I ) found that the outer edge of 



the gap (i.e. the inner edge of the circumbinary disk) depends on the binary eccentricity, 
with more eccentric binaries opening larger gaps. However, t he eccentricity itself increases 



Gold 


reich Sz Sari 


2003; 


Armitage & Nataraian 


2005; 


Cuadra et al. 


2009; 


Roedie et al. 


2011 


), resulting in a steady expansion of the gap. These two effect together would result 



in a limiting eccentricity of ~ 1, leading to a fast coalescence due to efficient gravitational 
wave emissions (see eq[2]). A smaller limiting eccentricity may result from a less efficient 
coupling between the gas and the binary, as expected if the disk moves farther out. Since 
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the evolution of e in an expanding gap has been performed with finite mass disks we 
postpone its discussion to the following. 

If the circumbinary disk is limited in mass, the evolution of the orbital separation and 
eccentricity can be quite different. Since in this scenario the disk is not continuously refilled 
from the outside, the gas mass initially available in the disk is the key parameter: 



• If the disk is S> M2, the evolution is similar to the infinite-mass disk, with the binary 
coalescing on a short time-scale. 

• If the mass in the circumbinary disk is less or of the order of M2, the interaction 
with the binary forces the whole circumbinary disk to move outward in few orbital 
periods. This expansion of the circumnuclear disk has been obs erved in simulations 
in which the compo nents of the binary have similar masses (e.g. 



Cuadra et al. 



2009 



Roedig et al. 



20111 ). In this case, the gas reaches distances > 4a, at which the 
interaction with the binary is not efficient anymore. A small amount of gas can then 
fall again closer to the binary because of orbital angular momentum exchange with the 
bulk of the disk, but most of the gas would never get close enough to significantly alter 
the evolution of the binary. Note that the expansion of the circumbinary disk is most 
effective for eccentric binaries. Binaries do not coalesce because of the interaction 
with too small M2), non- refuelled disks. 



In this simple description, the mass in gas within the disk is either accreted onto the 
MBHs or conserved. However, the gas in such a dense environment could be consumpted by 
star formation, decreasing the effective mass of the circumbinary disk. As a consequence, 
the disk could be initially ^ M2, but decreasing in mass with time, and could possibly fail 
in bringing the binary to the final coalescence. This scenario has been recently discussed in 



1 r 

Lodato et al. I (l2009l) . In this investigation stars are allowed to form in the disk whenever 
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it becomes gravitationally unstable. The rate of new star formation is obtained requiring 
that they would inject enough energy to keep the disk on the verge of fragmentation (i.e. 
providing an heating term exactly equal to the cooling losses in the disk). Even considering 
such a simple "thermal" feedback from the newly formed stars, the disk loses so much mass 
that the binary cannot reach the final coalescence, unless its initial separation is ^ 0.01 pc 
(for the MBH masses considered in the paper. Mi = 10^ Mq and M2 = 10^ M©). 

Other possible feedback terms, not included in this model, that can help in preventing 
such a strong gas consumption has been suggested by Lodato and collaborators, such 
as momentum feedback from stellar winds and supernovae explosions. Furthermore, the 
interaction between the MBHs and the forming gas clumps and stars, not considered in the 
investigation, could help the binary decay. The consumption of gas may not be a problem 
if large inflows of gas are present, as in the case of a continuously refuelled disk discussed 
above. 

In the finite r nass d isk scenario, the existence of a limiting eccentricity has been studied 



m 



Roedig et al. 



( 120111 ) through a suite of high resolution SPH simulations. They find a 
critical value Ccrit ~ 0.6 — 0.8. In these simulations, the initial ratio 6 between the gap size 
and the semimajor axis of the binary is 2, and can increase during the runs up to more 
than 4, when the interaction efficiency drops. The analytical model presented in the paper 
agrees with the simulations, predicting the limiting eccentricity to be: 



^crit 



0.66yin(5- 0.65) + 0.19. 



The initial choice of 5 = 2 is somewhat arbitrary. In reality, the feeding of a MBH 
binary forming in a gas-rich galaxy merger can be a very dynamic process, and the 
interaction with a single circumbinary disk could be too idealized a picture. Larger scale 
simulations show episodic gas i nflows due to the dynamica l evolution of the nucleus of the 



remnant (see e.g. 



Escala 



2006 



Hopkins fc Quataert 



2010). In this scenario the binary can 
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still interact with a disk and excavate a gap, but the size of it would be time- dependent 
(as in the simulations presented here), and would depend on the angular momentum 
distribution of the inflowing streams, resulting in a range of Ccrit- 

Note that the discussion above implicitly assumes that the MBHB and the 
circumnuclear disk corotate with each other. This is the natural outcome of a evolutionary 

le two MBHs orbit in a large scale 



2009al ). and open a gap in 



sequence in a gas rich galaxy mergers, in which t 
circumnuclear disk, are forced to corotate with it (jPotti et al. 
the very central region of the gas distribution. This picture, however, could not apply in 
gas-poorer mergers, or even in a gas rich scenario, if the circumnuclear disk formed during 
the merger fails in bringing the two MBHs to the final coalescence before it is consumed 
by star formation and/or MBH feedback. In one of these cases, an occasional small inflow 
of gas could happen with a random angular momentum, and could form a retrograde 
circumbinary disk, counterrotating with respect to the binary. 



Nixon et al. 



he evolution of a MBHB in a retrograde disk has been discussed in 
(I2OIII ). In this case, the gravitational interaction between the binary and the gas brakes 
both the components, so that, unlike in the prograde scenario, here the torques responsable 
for the binary shrinking and its eccentricity evolution causes the edge of the disk to move 
inwards. The binary, and the secondary MBH in particular, experiences the presence of a 
closer distribution of gas, which would imply a faster evolution, moving at a higher relative 
velocity, that results in a less effective interaction. Nixon and collaborators show that i) if 
the binary is initially not exactly circular (e < Ccrit ~ H/R, where H/R is the aspect ratio 
of the disk), and ii) if it interacts with the disk mainly at the apocenter, then the secondary 
evolves onto an almost radial orbit after interacting with a gas mass comparable to its own 



s 



■^Note that assuming the interaction to take place mainly at the apocenter is in agreement 
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The increase of the eccentricity to e ~ 1 in the retrograde case is due to direct accretion 
of hnear momentum from counterrotating material. Since the secondary has null radial 
velocity at the apocenter, before and after the interaction, and far from the apocenter the 
secondary is assumed to move on an unperturbed Keplerian orbit, the apocenter is constant. 
Interacting with counterrotating gas the secondary decreases its angular momentum, 
reducing its semimajor axis (of up to a factor of 2) and, most importantly, increasing its 
eccentricity (up to 1). At very high eccentricities the emission of gravitational waves can 
bring the binary to the final coalescence in less than an Hubble time (see eq. [2]). 

Note that this scenario suffers of the same disk consumption problem as the prograde 
one. If the disk is consumed by star formation or evacuated by MBH or supernovae 
feedbacks, the process stops. This makes this process particularly interesting for very 
unequal mass binaries (M2/M1 < H/R, less likely to form from galaxy mergers, as discussed 
in the Introduction). Fast inflows of gas, on timescales shorter than the consuption time, 
would help the coalescence in both the prograde and the retrograde scenario. 



4. Binary candidates in the realm of observations 



Despite being a natural outcome of galaxy mergers, MBH pairs are still elusive. 
Less than 20 systems with separations of ~ 10 pc to ~ 10 kpc pairs are of this kind are 
known to date. MBHs orbit in the common post-merger stellar environment, in-spiralling 
because of dynamical friction. They appear as a single galaxy (eventually, with disturbed 
morphology) with two active nuclei. Example s are the prototypica l case of NGC 6240 



(iKomossa et al. 



2OO3I ): the radio galaxy 3C75 (jHudson et al 



20061); the spiral galaxy NGC 



with the circumbinary disk to form after the MBHs bind in a close binary, as a consequence 
of a randomly oriented accretion event. 
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COSMOS J100043. 15+020637.2 fiComerford et al 



( I Green et al. 



3341 f lBarth et al. 1120081 ): the UL IRG Mrk463 fiBianchi et al. 1120081) : the interacting galaxy 



20091 ): and the quasar pair J1254+0846 



20101). All these objects have been discovered because of the presence of two 
resolved X-ray sources wandering in the merged galaxy. In order to look for these systems, 
an alternative app roach is to search for objects with tw o systems of narrow lines at slightly 



different redshifts (IWang et al. 



2009 



the Sloan Digital Sky Survey (SDSS 



Liu et al. 



York et al. 



2010 ). Large spectroscopic surveys, like 



2000l ). have been used to search for these 



systems. Follow-up observations were then used to discriminate be tween dual AGN an d 



single AGN wi th complex gas dynamics in the narrow line region (jWang et al. 



2009; 



Fu et al. 



20111 ). 



At separations of < 10 pc the two MBHs start experiencing their own gravitational 
interaction, binding in a binary. These systems cannot be spatially resolved in optical 
and X-ray observations, and radio i nterferometry is required. This has been su ccessfully 



done only in the case of 0402-1-379 (IManess et al. 



2004 



Rodriguez et al. 



20061 ). The two 



fiat-spectrum radio sources, corresponding to the two components of the candidate MBHB, 
have a projected separation of ~ 7 pc (few milliarcsec aX z = 0.055). We note however that 
radio interferometry at very high spatial resolution is not an efficient technique to search 
for rare objects as MBHBs, because of the li mited field of views a nd the requirement that 



both the MBHs are radio- luminous (see, f.i.. 



Burke-Spolaor 



201 ih 



Another approach to look for MBHBs is to study periodic variations in the luminosity 
of some AGN. Th e only MBHB candida te selected on these bases up to now is the BL Lac 



object OJ287 (see 



Valtonen et. al 



2008 



and references therein). It shows a 11 yr periodic 
flaring. In the MBHB scenario, the secondary MBH orbits on a tilted plane with respect to 
the accrtion disk of the primary. The flares are associated to the passages of the secondary 
MBH through the nodes. However, this is not the only available explanation of the peculiar 
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behavior of this source (jVillforth et ah Il2010l ). 



More promisingly, signatures of MBHBs have been searched for in optical and NIR 
spectroscopic databases. According to the MBHB hypothesis, if at least one of the MBHs 
is active, the broad lines (BLs) emitted by gas bound to it may be red- or blue-shifted 
with respect to their host galaxy sy stemic recessional velocity, as a con sequence of 



the Keplerian motion of the binary ( iBegelman. Blandford &: Rees 



1980). Therefore, 



looking for quasars with significant velocity shifts (>few hundreds km s~^) can be a valid 
approach to systematically search for MBHB candidates over large fields. This technique 
does not suffer of any angular resolution limitations: Actually, the closer (and more 
massive) the binary is, the more shifted/deformed the BLs are. Five MBHB candidat es 



have 



2008 



jeen individually found in this way: SPSS J0927-I-2943 



Bogdanovic. 



Boroson fc Lauer 



(Decarli et al. 



Eracleous &: Sigurdsson 



200S 



Dotti et al. 



Komossa. Zhou fc Lu 
2009bi. J153 6+0441 



20091), J1050+3456 fIShields. Bonning fc Salviander 



2010h . and J0932-F0318 flBarrows et al 



201 ih 



20091). JlOOO +2233 



Tsalmantza et al. 



(120111 ) recently applied this technique in a systematic way over the whole spectroscopic 



sample of SPSS, using a method developed for searching unresolved gravitational 



lenses (ITsalmantza fc Hogg 



20111 ). This analysis resulted in 4 new MBHB candidates 

1714+3327) and the co nfirmation of all the 



Eracleous et al. 



2011 



searched fo r objects 



2010) ■ Among 



(J1012+2613, J1154+0134, J1539+3333 and 
previously known objects. In a similar study, 

with anomalous H/3 profiles in the SPSS quasar catalogue (iSchneider et al 
them, they identified 88 sources showing velocity shifts between the broad H/3 line peak 
and the rest-frame of the narrow emission lines. 

It should be noted that alternative interpretations for the spectral properties of the 
known candidates are available: 



i- Modest line shifts f< 500 km s ^) are often observed in "normal" AGN 
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Velocity [xlOOO km s"'] 



Fig. 3. — Velocity diagrams of Ha (red, solid lines), H/3 (blue, dotted lines) and Mgll (green, 
dashed lines) for all the spectroscopically-identified MB HB candidates. The flux o f H/3 is 
scaled up to match the one of Ha or Mgll. Adapted from iTsalmantza et al. I ( 120111 ). 
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(iBonning. Shields &: Salviander 1120071 ). 



ii- Similarly, small velocity shifts (<4000 km s can be associated to the remnant 
of a binary coalescence, reco iling because of anisotropic gravitational wave emission 



( iKomossa. Zhou fc Lu 



2008 



in- An unobscured MBHB with both M BHs active could resamble t he spectrum 
of a double peaked emitter (see, e.g.. 



Eracleous Sz Halpern 



19941 ). where broad 



doublepeaked lines are emitted because of the almost edgeon, disklike structure of the 
broad line region of a single MBH. 

iv- Chance superposition of two AGN (or an AGN-galaxy superposition) within the 



angular resol ution of the usee 



line systems ( iHeckman et al. 



spect rograph can also mimic velocity shifts of different 



20091). 



The simplest way to discriminate between these scenarios and the MBHB hypothesis 
would be to look for a periodic oscillation of the broad line shifts around the host galaxy 



redshift. However, the orbital 



jeriod of the binary could be too long to be easily observed 



( jBegelman. Blandford fc Rees 



Eracleous et al. 



(I2OIII ) observed a 



19801 1. Noticeably, 

variation in the shifts at two different epochs in 14 out of 88 candidates, with resulting 
accelerations between -120 and +120 km/s/yr. Longer temporal baselines are needed to 
prove the MBH binary interpretation for these objects. 

In order to increase the number of known MBHB candidates, and to confirm their 
interpretation, it is therefore of fundamental importance to identify new signatures of 



^However, if t 



le galaxy merger is gas rich, the maximum recoil 



2010 



be < 100 km s^ (IBogdanovic. Reynolds fc Miller 



Kesden. Sperhake fc Berti 



2OI0I ) 



2007 



veloci t y is expected to 



Dotti et al. 



2010 



Volonteri et al. 
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MBHBs. The simultaneous observation of various MBHB signatures could represent the 
only way to firmly validate the MBHB scenario in the known candidates. A possibility is 
to look for peculiar flux ratios between high- and low-ionization broad lines. The broad 
line region of each MBH in a binary can be perturbed and disrupted by the gravitational 
potential of the companion. External shells of the broad line region (where most of 
the low-ionization line flux is emitted) are affected first, resulting in peculiar flux ratios 



( iMontuori et al. 



20111 ). This criterion is particularly interesting for quasars at high redshift 



(z ~ 2), where high and low ionization lines are observable in large surveys as the SDSS, 
while the most prominent narrow lines, needed to measure a shift related to the orbital 
motions, are not present in the SDSS spectra anymore. 

At even closer separations between the two MBHs, when the size of the BLR is 
significantly larger than the semimajor axis of the binary (a < 0.01 pc), the optical and 
UV spectral features disc ussed above become more complex and not directly related with 



the period of the binary (IShen &: Loeb 



2010l ). However, in this case typical MBHB periods 



are ^ 10 yr, opening the interesting opportunity of directly det ect periodic variability 
of th e system, related to the periodicity of the accretion fl ows (lArtymowicz fc Lubow 



1996 



Hayasaki. Mineshige fc Ho 



2008 



Sesana et al. 



20111 ). Such close separations are 



particularly in teresting, since they will be proven b y the ongoing and fu ture pulsar timing 



arrays (PTAs, 



Hobbs et al. 



2010l ). In this context. 



Sesana et al. 



(120 111 ) estimate that up 



to few hundred MBHBs contributing to the GW signal in the PTA band may be identified 
through their periodicity in future X-ray all sky surveys. Among those, few exceptionally 
bright sources may be resolved both in the GW and in the electromagnetic window through 
the detection of peculiar double K-a fluorescence lines, offering a unique multimessenger 



astron omy opportunity. An alternative possibihty, suggested by 



Tanaka. Menou fc Haiman 



( l201ll ). is that the presence of a circumbinary cavity results in a suppression of the UV 



soft-X emission. MBHBs close to coalescence may therefore be identified as AGNs with 
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exceptionally faint UV X-ray continuum. 



Conclusions 



We reviewed the most recent findings about the dynamical evolution of MBHBs and 
their detectability. Regarding the binary dynamics, in the last few years we recognized 
the importance of the medium/large scale galactic structures 100 pc) in the dynamical 
evolution of the binaries. In gas free environments, the shape of the bulge potential is 
directly related to the possibility of a MBHB to reach the final coalescence. In a spherical 
system, the stars that can interact with the binary (i.e. that orb it within its loss-cone) 



Milosavljevic fc Merritt 



are e vacuated from the center before the binary can coalesce (e.g. 
20011 ). A fast refilling of the loss-cone, that can result in the merger of the two MBHs, is 
possible in triaxial systems, in which the angular momentum components of the orbiting 



stars are not conserved (e.g. 



Poon&Merritt 



20041 ). Recently, thanks to the advent of GPU 



computing, large s cale galaxy mergers proved the occurrence o : 



realistic scenarios (IKhan. Just fc Merritt 



2011 



Preto et al. 



such a replenishment in mor e 



2011 



Gualandris fc Merritt 



20111 ). The evolution of the binary eccentricity e depends on the dynamical properties of the 
core as well. If the MBHs are embedded in a non rota ting stellar system, the general trend 



is towards an increase of e with time (e.g. 



Quinlan 



19961 ). This promotes the coalescence 



of the MBHs, since gravitational wave emission is more efficient for eccentric binaries. In 
rotating systems, on the other hand, the evolution of e depends on the orientation of the 



binary: a binary co-rotating with the stellar cusp t ends to decrease its eccentricity, w 



lile 



201l|). 



in the counter-rotating case e grows up to > 0.95 (iSesana. Gualandris fc Dotti 
This clear cut scenario could be modified by the interaction with stars on quasi-radial, 
centrophilic orbits. A study of the orbital properties of these stars has not been presented 
to date. 
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Galaxy mergers can easily promote strong inflows of gas towards the center of the 



galaxy remnant (e.g. 



Mayer et al. 



20071 ). Hence, it is fundamental to understand how the 
presence of massive gas structures in the cores modify the dynamical evolution of the 
forming MBHBs. It has been proven that the interaction with circumnuclear (~ 1 00 pc) 



disks can result in a fast (< 10^ yr) formation of a MBHB (e.g. 



Escala et al. 



20051). After 



this fast transient, the bi nary is thought to open a gap in the central gas distribution (e.g. 



Lin &: Papaloizou 



1979al ). and any further evolution is mediated by the interaction between 



the MBHs and the inner edge of the circumbinary disk. 

Simulations and analytical studies about the interaction of MBHBs and circumbinary 
disks have improved our knowledge of the physical processes in play, and their effect onto 
the binary evolution. Howerver, a complete understanding of the binary/disk interaction is 
still to come. For example, in many investigations the circumbinary disk is assumed to be 
cylindrically symmetric, i.e., the study is reduced to an effective one-dimensional problem. 
The assumption of cylindrical symmetry removes any possibility of studying the effects of 
structures in the disk, and, most importantly, of gas streams p eriodically inflowing from the 



disk onto the two MBHs, routinely 



Havasaki. Mineshige fc Ho 



2008 



observed in simulations (e.g. 



Cuadra et al. 



2009 



Artvrnowicz fc Lubow 



Roedig et al. 



2011 



1996 



Sesana et al. 



201l|). The torques exerted by th ese inflowing streams has not been studied in detail 



yet (with the notable exception of 



MacFadyen Sz Milosavljevic 



20081), and could provide 



additional help (or resistance) in bringing the binary to the final coalescence. 

The final fate of a binary embedded in such a circumbinary disk is still debated. If the 
disk is continously refueled from any larger-scale gas distribution, a fast coalescence can 
easely be achieved. However, if the binary cannot interact with enough gas (e.g. because 



it turns in stars), t he circumbinary dis 



to coalescence (e.g. 



Lodato et al. 



i gets evacuated and fails in bringing the MBHs 



20091 ). As a consequence, as in the stellar scenario, the 
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final fate of the binary depends on the properties of larger scale structures, and its abitily 
to efficiently refuel the proximity of the MBHB with fresh gas. In principle, in presence 
of an intense inflow toward the center, the binary could fail in opening a gap at all, and 
the interaction between M BHs and a closer /d enser gas structure could result in a faster 



coalescence of the MBHB (lEscala et al. 



20051 ). Furthermore, if the angular momentum of 



gas can be efficiently re-shuffied, inflowing streams co uld form counter-rot ating circumbinary 



disks, that can also promote the binary coalescence (INixon et al 



the formation of a gap has been observed in large scale simulations (lEscala et al. 



201 If). Only recently 



2005 



Dotti et al. 



2009al ). in which the evolution of an extended (~ 100 pc) massive disk is 
followed (as massive as the binary in Escala et al., up to ^ IOMbhb in Dotti et al., see 
flgure|2]). The spatial and mass resolution of these simulations do not allow yet a detailed 
study of the sub-pc evolution of the binary, down to a possible coalescence. Simulating 
the evolution of a sub-pc binary starting from large scale initial conditions, that can 
constrain the properties of the nuclear inflows together with the evolution of the binary, 
is the fundamental improvement needed to build a coherent picture of MBHBs in gas rich 
environment. 

To summarize the recent flndings present in literature, we can draw a comparison 
between orbital decay timescales obtained considering different scenarios, for equal mass 
binaries. 



In dense stellar environments: z/the loss-cone of the MBHB is constantly refllled (see 
Section 2.1), a binary of lO^M© with an initial separation of ao ~ 1 pc will coalesce in 
~ 3 X 10'' — 10^ years , while a binary 100 times more r nassive ^yill inspiral for abou t 



10^ — 10^ years before reaching the flnal coalescence (jSesana 



2010; 



Preto et al. 



201l|). 



in gas rich environments: if the interaction with a steady, long-lived, corotating, 
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maximally massive circumbinary disk is responsable for the MBHB orbital decay (see 
Section 3.2), the orbital decay timescale for a IO^Mq equal mass binary is less than 
the age of the Universe if the binary starts at an initial separation ag ~ 0.03 pc (~ 10 
times smaller than the separation at which the MBHs bind in a bina ry aBHB), and is 



2 orders of magnitude shorter for ao ~ 0.01 pc ( jCuadra et al 



20091 ). For a lO^M© 



equal mass binary, t 



( ICuadra et al. 



20091 ) 



lis timescale is less than the age of the Universe if ~ 0.05 pc 
100 times less than aeHB- 



• in gas rich environments: if the MBHB interacts with a continuous sequence of 
counterrotating accretion disks with an accretion rate corresponding to the MBH 
Eddington limit, the orbital decay timescale is of the order of 10*^ yr, regardless the 
binary mass, for < ctBHB- In this case, the coalescence timescale increases linearly 
with the inverse of the accretion rate. 



Note that the estimates in gas rich environment should be considered as lower limits, 
since they assume continuous accretion and, as stressed above, a continuous re-fuelling of 
the disks from larger scales. A single, not re-fuelled disk with enough mass and angular 
momentum to bring the MBHs to the final coale scence form an ^ 0.1 pc would undergo 



fragmentation and star formation, as discussed in 



Lodato et al. 



(120091 ) 



The presence of gas close to the binary is necessary to its detection. If at least one of the 
two MBHs is active, the orbital velocity of the binary can res ult in a frequency shift between 



the b road emission lines and the narrow emission lines (e.g. 



Begelman. Blandford fc Rees 



1980l ). This shift is expected to change periodically, on the orbi tal period. Using this 



technique few tens of MBHB candidates have been selected (e.g. 



Eracleous et al. 



Tsalmantza et al. 



2011 



201ll ). For all the candidates discussed to date, possible explanations other 



than a MBHB have been proposed (see section Hj). Moreover, the orbital period expected 
for such binaries is often too long to be observed, thus the periodic variation of the velocity 
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shift cannot be used (yet) to determine the real nature of the candidates. It is therefore 
fundamental to couple this with other (independent) signatures of MBHBs, to confirm their 
nature. 

In the near future, space based interferometers like NGO will detect the GWs emitted 
in the late inspiral and final coalescence of MBHBs, opening the fascinating prospective 
of multimessenger astronomy. The identification of an electromagnetic counterpart to a 
gravitational wave detection of merging MBHs can teach us about relativity, accretion 



physics, galaxy formation and evolution and, more in ge neral, cosmo. 



of potential EM signatures have been proposed (see, e.g. 



Schnittman 



ogy. A large variety 



2011 



and references 



therein). In particular, thanks to the recent quick progresses in numerical relativity. 



simulations o 



Bode et al. 



the coalescence process in presence of matter is now becoming possib 



2010 



Moesta et al. 



Palenzuela. Lehner &: Liebling 



2010 



Farris. Liu &: Shapiro 



2011 



20111 ). and in the coming years may provide valuable insights about the 
nature of the associated electromagnetic signals. However, most of the counterparts 
discussed to date depend on some assumptions on the distribution of gas and/or stars in the 
immediate vicinities of the MBHB. A theoretical comprehension of the dynamics of matter 
close to the binary, necessary to understand the fate of a MBHB, is also fundamental to 
constrain the observability of an electromagnetic counterpart. 

Tracing the formation, evolution and final fate of MBHBs is certainly one of the open 
challenges of contemporary astrophysics. A better understanding of the interaction between 
MBHs and gas, and the prediction of new peculiar observational features of MBHBs are 
needed to unambiguosly constrain their properties and demography, adding an important 
missing piece to the galaxy evolution puzzle. 



We are grateful to Jorge Cuadra, Monica Colpi and Paraskevi Tsalmantza for useful 
discussions and comments. 
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